X-RAY EMISSION OF YOUNG SN lA REMNANTS 
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We present results of hydrodynamical simulations of young supernova remnants. To model the ejecta, we 
J> ■ use several models (discussed in literature) of type la supernova explosions with different abundances. Our 
0^ . hydro models are one-dimensional and spherically symmetrical, but they take into account ionization kinetics 
I with all important processes. We include detailed calculations for the X-ray emission, allowing for time- 
CsJ dependent ionization and recombination. In particular, we compare the computed X-ray spectra with recent 
. XMM-Newton observations of the Tycho SN remnant. Our goal is to find the most viable thermonuclear 
Q \ SN model that gives good fits to both these X-ray observations and typical SN la light curves. 

^: 

O ' Introduction 

Recent observations with XMM-Newton telescope with high spectral and space resolution, give us an 

■ opportunity to get clearer understanding of nature and origin of X-ray sources. For example, with new 
J> \ data on supernova remnants we can find out more about the supernova explosion itself. So far, there exist 

many explosion models proposed by theorists for different types of supernovae, but still there are no definite 
criteria to decide which of the models are realized in nature. During the first months after an explosion one 
5^ \ can examine a theoretical model by calculating bolometric and monochromatic light curves and spectra. 
Later on, gas in the ejecta cools down and becomes almost unobservable. The next opportunity to analyze 
the ejecta is on the stage of a young supernova remnant (SNR), when noticeable amount of circumstellar gas 
is swept up. At this moment a reverse shock forms, goes inwards the ejecta and illuminates it once again. 

We will focus only on Type la Supernovae (SN la) and compare the results for four theoretical models of 
explosion. We are eager to find a model which would be in good agreement with observations of SN la both 
at the epoch of maximum light and on the young remnant stages. 

In this work, we calculate hydrodynamical evolution and X-ray emission of a supernova remnant at the 
age of 430 years, which corresponds to the age of Tycho SNR, and compare the results with spectra and 
images of this remnant obtained by XMM-Newton space telescope (Decourchelle et al. 2001). We examine 
the same models that where applied to our broad-band light curve modelling (Blinnikov, Sorokina 2003), so 
we will be able to judge on the correctness of these models from two points of view. 

Quite a while ago, it was concluded that Tycho's supernova was of Type I (Baade 1945). After that time 
there were many other papers which suggested other types for Tycho. But the morphology of the remnant, 
its shell-like structure, says that most probably it is still the remnant of Type la supernova. 

The basic parameters of the Tycho remnant are the following: the age is 430 years, the angular size is 
~ 8', thus the radius is ~ 2 — 3 pc (depending on adopted distance). It has almost spherical shape, so we 
do not need any additional assumptions to model it with a ID code. 

Tycho's remnant was observed for many years in radio wavelength and optics with terrestrial instruments 
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(e.g., van den Bergh 1971, van den Bergh et al. 1973, Duin, Strom 1975, Dickel et al. 1982) and in X- 
rays with different space telescopes (e.g., Davison et al. 1976, Seward et al. 1983, Tsunemi et al. 1986, 
Hwang, Gotthelf 1997, Decourchelle et al. 2001). Its images and spectra were analyzed in many papers 
(e.g., Chevalier, Raymond 1978, Hamilton et al. 1986, Itoh et al. 1988, Brinkmann et al. 1989, Vancura et 
al. 1995, Hwang et al. 1998). The instruments become better every decade. The most recent observations 
with XMM-Newton provide us with data of excellent spatial resolution, which allows us to examine not only 
crude models to reproduce the emission of the remnant, but also to change between models which differ not 
very strongly. 

Models 

To simulate supernova ejecta in our calculations of evolution and X-ray emission of remnant we have 
chosen four SN la explosion models. Three of them are Chandrasekhar-mass models and one is a sub- 
Chandrasekhar one: 

• the classical deflagration model W7 (Nomoto et al. 1984). Mwd = 1.38M0, Eq = 1.2 x 10^^ ergs, 
M(56Ni) = O.6M0. 

• the deflagration model by Reinecke et al. (2002) (hereafter MRO). Mwd = 1.38M0, Eq = 4.6 x 

10^° ergs, M(5*'Ni) = O.43M0. 

• the delayed detonation model DD4 (Woosley & Weaver 1994b). Mwd = 1.38M0, Eq = 1.2 x 10^^ ergs, 
M(56Ni) = O.6M0. 

• the sub- Chandrasekhar-mass detonation model with low ^^Ni production (hereafter, WD065; Ruiz- 
Lapuente et al. 1993). Mwd = O.65M0, Eo = 5.6 x 10^° ergs, M{^^m) = 0.5Mq. 

W7, DD4, and WD065 are well known SN la models. The last one is the best for modelling light curves 
and late spectra of very sub-luminous SNe la, like 1991bg. Two others are usually considered as the best 
models for typical SN la. 

The recent MRO model is originally 3D, but we have averaged it over the whole Air, since our hydro code 
is ID. MRO is almost a "first principle" model, and uses much less free parameters for flame modelling than 
all previous calculations. It is less energetic, contains less ^^Ni, but the latter is located in the outermost 
layers of ejecta due to 3D Rayleigh-Taylor instability. The combination of all these factors leads to the 
situation when the outermost ^^Ni layer has similar velocities for all Chandrasekhar-mass models. 

Prom our light curve modelling (Blinnikov, Sorokina 2003) we have found that MRO model is the best 
for UBVI bands, though the bolometric light curve is more similar to observations for W7. We discussed 
there that a bit more energetic model, and equally mixed as MRO, would be the best from the light curve 
modelling point of view. 

We surround ejecta by a motionless gas of constant density (5 x 10~^^ g/cm^), constant temperature 
(10^ K), and with Solar abundances. We start calculations of the remnant evolution at the age of ejecta 
from several years to several decades. 

Method 

To model the hydrodynamical evolution of the supernova remnant we use hydro code from the pack- 
age STELLA (Blinnikov et al. 2000). To find out the X-ray spectra we use the code written by P.L. for 
calculations of collisional ionization of a stationary plasma, which takes into account basic processes: ioniza- 
tion by electron impact, autoionization, photorecombination, dielcctronic recombination, charge transfer. It 
calculates bremsstrahlung (free- free), free-bound, bound-bound, two-photon emission. This code has been 
elaborated by S.B. and E.S. into a time-dependent variant (n^ = f{T(t),n(t),Z,...) for all ionic species). 
Input parameters are current temperature, density, abundance. 

The models were calculated as follows. We simulate an SN expansion into CSM, assuming that the remnant 
is transparent, i.e. radiation does not affect remnant's dynamics, so we can calculate hydrodynamics and 
radiation separately. During the calculation we record history of temperature and density variations for 
each mesh zone and then we use this data to calculate time-dependent ionization. The resulting ionization 
stages for all elements at the age of Tycho remnant are used to evaluate an X-ray spectrum. 

The problem here can arise in the correct hydro calculations without knowing the ionization stage of gas. 
For CSM, if it consists mostly of hydrogen, ionization should not affect hydro evolution strongly. But the 
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ejecta is more metal- abundant. Number of free electrons can vary by a factor of a few for different ionization 
stages, that leads to strong changes in pressure within the ejecta and, therefore, to possible differences in 
the speed of the reverse shock. Because of limited computational resources we decide to start with flow 
models with a simplified equation of state. We assume a fixed, but not very low or very high, ionization 
stage (say, the 8th) of shocked gas during the calculations. In future we plan to iterate the ionization stage 
and repeat hydro calculations, which should make them more self-consistent. 

Results and comparison with observations of Tycho remnant 

The resulting hydro dynamical structure of our models at the age of Tycho remnant (430 years) is shown 
in Figjl] and Fig.|2|. In each model the reverse shock is formed and goes inward the ejecta, but does not reach 
the center yet. 
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Fig. 1. Temperature (solid), density (dash-dotted) and velocity (dashed) profiles for models: left: WD065; 
right: MRO 
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Fig. 2. Temperature (solid), density (dash-dotted) and velocity (dashed) profiles for models: left: DD4; right: 
W7 

Due to different explosion energy in our models, they expand to the given age till different radii, the 
position of the contact discontinuity (that separates ejecta and ambient medium) for different explosion 
models also varies from 1.8 to 2.5 pc. Maximum value of the shocked ejecta density can vary in a range of 
an order of magnitude for different models. The differences in temperature of heated ejecta are of the order 
of half magnitude for different models. The width of the heated (therefore, emitting) region varies also quite 
noticeably. All these varieties lead to the differences in the emission one could observe from the remnant. 



4 



Results of the spectra simulations for different explosion models are shown in Fig.^ and Fig.^. The gray 
line corresponds to Nh = 4 x 10^^ cm~^ galaxy column density, the black one - to Nu = 9 x 10^^ cm~^. 

These plots show that we can rule out the WD065 model (it shows very weak Fe K-line, contrary to 
the EPIC observation with XMM-Newton by Decourchelle et al. 2001). All Chandrasekhar-mass models 
produce almost the same spectra with similar strength of lines. Though Fe K emission is a bit stronger in 
MRO, it is not possible to judge which of these models is better. 
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Fig. 3. Full theoretical X-ray spectra (with the account of interstellar absorption for column density Nh) of the 
simulated remnant in XMM-Newton range (0.2-10 keV) for various types of models: left: WD065; right: MRO 
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Fig. 4. Full theoretical X-ray spectra (with the account of interstellar absorption for column density Nh) of the 
simulated remnant in XMM-Newton range (0.2-10 keV) for various types of models: left: DD4; right: W7. 



Much more drastic differences arise when we look at the profiles of the remnant in bandwidths of various 
ion lines and compare these with observed azimuthally averaged radial profiles of the deconvolved images 
in various lines (Decourchelle et al. 2001). 

We plot (Fig.|5| and Fig.^) the profiles of the simulated remnant in different lines (e.g., in FeXVII line 
and Si K line). Models DD4 and W7 give comparable profiles, but Si and Fe lines are separated strongly in 
space, WD065 is distinguished strongly from the others and from the observational data, while MRO shows 
the best agreement with observations: the diameters of the remnant in Si and Fe lines are very close, an 
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amplitude of the luminosity variations in the corresponding lines is also more similar to what is observed by 
Decourchelle et al. 2001. 

Model: p800_5d24; Age: 430 yr Model: mr_800_5d24; Age: 430 yr 
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Fig. 5. Profiles of the remnant's brightness in lines of Si K and FeXVII ions, integrated over line of sight. Dashed 
line shows a profile of emission within a range 1.67-2.0 keV (Si K) and 775-855 eV (FeXVII). Solid line shows 
emission within a range of one bin {AE = 13.6 eV 1 keV) that presumably corresponds to the emission line of 
Si K or FeXVII. Left: WD065; right: MRO 
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Fig. 6. Profiles of the remnant's brightness in lines of Si K and FeXVII ions, integrated over line of sight. Dashed 
line shows a profile of emission within a range 1.67-2.0 keV (Si K) and 775-855 eV (FeXVII). Solid line shows 
emission within a range of one bin (AE = 13.6 eV 1 keV) that presumably corresponds to the emission line of 
Si K or FeXVII. left: DD4; right: W7 



Discussion and prospects 

Though the results are very preliminary, we can rule out the sub-Chandrasekhar-mass model WD065 for 
the Tycho SNR. The most popular models DD4 and W7 show relatively good agreement with observational 
data (though the radial profiles are not very encouraging). Nevertheless, the most promising model is still 
MRO. It fits better both the observed X-ray spectrum and the brightness profiles in different lines. 

To make our results more reliable, we still need to upgrade our techniques in order to obtain a better 
agreement with X-ray observations. There is a number of issues that have to be worked out. 
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Since we do not know exactly the ionization stage (thus the "true" equation of state) of plasma on the 
hydro calculation step, we should calculate it by our time-dependent code. But due to our limited computer 
resources we cannot make it in a fully self-consistent way immediately. Instead, we can calculate separately 
and save the history of ionization and then perform the second iteration of flow calculation with the new 
(updated) set of ionization data. 

In the modelling we assume Tg = Tj in the calculations. We can introduce a parameter that specifies 
relationship between Tg and Tj. Variations of this parameter will affect the spectrum shape. Thus we can 
find a value of the parameter that will make a best fit to spectrum. 

Wc have only tested CSM with a constant density. In near future we plan to test also "windy" CSM, with 
density gradient. That will change the dynamical evolution and, probably, the emission of our models. 

At last, looking at the spectra one can see, that Decourchelle et al. 2001 show very broad lines and we have 
them rather narrow, although we have taken into account thermal and Doppler broadening. Convolving our 
theoretical curves with XMM-Newton response matrix (EPIC, MOS and PN instruments) would give us a 
possibility for more accurate data fitting. 

Taking into account all the above issues we will make the comparison of models and observations more 
correct and fruitful. Nevertheless, we believe that the main conclusion of this work will remain unchanged: 
the new 3D SN la model MRO, due to its low energetics and a strong mixing, seems the best to model Sn la 
remnants, as well as their light curves. 
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